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Abstract. We give an entirely new approach to the problem of 
mutually unbiased bases (MUBs), based on a Fourier analytic tech- 
nique in additive combinatorics. The method provides a short and 
elegant generalization of the fact that there are at most d+ 1 MUBs 
in C d . It may also yield a proof that no complete system of MUBs 
exists in some composite dimensions - a long standing open prob- 
lem. 
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1. Introduction 

In this paper we introduce a novel approach to the problem of mu- 
tually unbiased bases in C d . Surprisingly enough, the required Fourier 
analytic technique is borrowed from additive combinatorics - a seem- 
ingly unrelated branch of mathematics. 

The paper is organized as follows. In the Introduction we recall some 
basic notions and results concerning mutually unbiased bases (MUBs). 
In Section [2] we describe how the problem of MUBs fits into a general 
scheme in additive combinatorics - a scheme we will call Delsarte's 
method. We then apply this method to prove Theorem I2.2[ an elegant 
generalization of the fact that there are at most d + 1 MUBs in C d . 
Finally, in Section [3] we indicate the limitations of the method by intro- 
ducing the notion of pseudo-MUBs, and discuss the possible existence 
of such in the case d = 6. 

Recall that given an orthonormal basis A = {e 1; . . . , e^} in C d , a 

unit vector v is called unbiased to A if |(v, e k )\ = —= for all 1 < 

v d 
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k < d. Two orthonormal bases in C d , A = {ei,...,ed} and B = 
{fi, . . . , id} are called unbiased if for every 1 < j, k < d, \ (e,, i k ) \ = —j=. 

A collection B ,...B m of orthonormal bases is said to be (pairwise) 
mutually unbiased if every two of them are unbiased. What is the 
maximal number of pairwise mutually unbiased bases (MUBs) in C d ? 
This question originates from quantum information theory and has 
been investigated thoroughly over the past decades (see [15] for a recent 
comprehensive survey on MUBs). The following result is well-known 



Theorem 1.1. The number of mutually unbiased bases in C d cannot 
exceed d + 1 . 

We will generalize this fact in Theorem 12.21 below. The other im- 
portant well-known result concerns prime-power dimensions (see e.g. 



Theorem 1.2. A collection ofd + 1 mutually unbiased bases (called a 
complete set of MUBs) can be constructed if the dimension d is a prime 
or a prime-power. 

However, if the dimension d = p" 1 . . . p^ k is composite then very little 
is known except for the fact that there are at least p- J + 1 mutually 
unbiased bases in C d where p™ 3 is the smallest of the prime-power 
divisors. In some specific square dimensions there is also a construction 
based on orthogonal Latin squares which yields more MUBs thanp" J +l 
(see |26j). The following basic problem, however, remains open for all 
non-primepower dimensions: 

Problem 1.3. Does a complete set ofd+1 mutually unbiased bases 
exist in C d if d is not a prime-power? 

The answer is not known even for d = 6, despite considerable efforts 
over the past few years ([U El CH ESI ED E2 E])- The case d = 6 is 
particularly tempting because it seems to be the simplest to handle with 
algebraic and numerical methods. As of now, some infinite families of 
MUB-triplets in C 6 have been constructed ([281 HH]), but numerical 
evidence suggests that there exist no MUB-quartets [SI [71 [91 [28] . 

It will also be important for us to recall that mutually unbiased 
bases are naturally related to complex Hadamard matrices. Indeed, 
if the bases Bo, ... , B m are mutually unbiased we may identify each 
B\ = {e^, . . . , e^} with the unitary matrix 



(see e.g. HJ H [27]): 



[H HH H21 HS1 HZl [201 [27]). 




A FOURIER ANALYTIC APPROACH TO MUBS 



3 



i. e. the k-th column of Hi consists of the coordinates of the k-th vector 
of E\ in the basis B . (Throughout the paper the scalar product (., .) 
of is conjugate-linear in the first variable and linear in the second.) 
With this convention, H Q = I the identity matrix and all other matrices 
are unitary and have entries of modulus l/y/d. Therefore, the matrices 
H[ = y/dHi have all entries of modulus 1 and complex orthogonal rows 
(and columns). Such matrices are called complex Hadamard matrices. 
It is thus clear that the existence of a family of mutually unbiased 
bases Bo, ... , B m is equivalent to the existence of a family of complex 
Hadamard matrices H' 1 ,...,H' m such that for all 1 < j ^ k < m, 
-^qH*H' k is again a complex Hadamard matrix. In such a case we will 
say that these complex Hadamard matrices are mutually unbiased. 

A complete classification of MUBs up to dimension 5 (see [5]) is 
based on the classification of complex Hadamard matrices (see [IS])- 
However, the classification of complex Hadamard matrices in dimension 
6 is still out of reach despite recent efforts [21 I2U [211251 CE2] • 

In this paper we will use the above connection of MUBs to complex 
Hadamard matrices to apply a Fourier analytic approach, borrowed 
from additive combinatorics. 

2. Mutually unbiased bases, difference sets and 
Dels arte' s method 

In this section we describe a general scheme in additive combina- 
torics, and show how the problem of mutually unbiased bases fit into 
this scheme. 

Let G be a compact Abelian group, and let a symmetric subset A = 
—A C G, G A be given. We will call A the 'forbidden' set. We would 
like to determine the maximal cardinality of a set B = {b\, . . . b m } C G 
such that all differences bj—bk € A c U{0} (in other words, all differences 
avoid the forbidden set A). Some well-known examples of this general 
scheme are present in coding theory ([II]), sphere-packings (P2]), and 
sets avoiding square differences in number theory ([22]). 

We now describe a general method to tackle such problems. To the 
best of my knowledge it was first introduced by Delsarte (in a less 
general terminology) in connection with binary codes with prescribed 
Hamming distance. The method is also 'folklore' in the additive com- 
binatorics community and I was introduced to it by Imre Z. Ruzsa 

([21). 
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We are looking for a 'witness' function h : G — > M with the following 
properties. 

• his an even function, h(x) = h(—x), such that the Fourier inversion 
formula holds for h (in particular, h can be any finite linear combination 
of characters on G). 

• h(x) < for all x G A c 

• Kl) > for all 7 G G 

• h(0) = 1. 

Lemma 2.1. (Delsarte's method) 

Given a function h : G — > K. wnt/i £/ie properties above, we can conclude 
that for any B = {6 ls . . .6 m } C G snc/i i/iaf bj — b k E A c U {0} t/ie 
cardinality of B is bounded by \B\ < h(0). 

Proof. For any 7 G G define -8(7) = ^2T=ij(bj)- Now, evaluate 

(1) s = J2\b(i)\ 2 Hi)- 

All terms are nonnegative, and the term corresponding to 7 = (the 
trivial character, i.e. 7(2;) = 1 for all x G G) gives \B(0)\ 2 h(0) = \B\ 2 . 
Therefore 

(2) S>\B\ 2 . 

On the other hand, (5(7) | 2 = X/jfc7(fy — ^fc)' anc ^ therefore 5 = 
Y^jjk^i^i ~ bk)h(j). Summing up for fixed j, k we get 
S 7 7(frj — bk)h(y) = h(bj — bk) (the Fourier inversion formula), and 
therefore S = ^2j k h(bj — bk). Notice that j = k happens |5|-many 
times, and all the other terms (when j 7^ k) are non-positive because 
bj — bk G A c , and h is required to be non-positive there. Therefore 

(3) S< h(0)\B\. 

Comparing the two estimates (|2J), © we obtain \B\ < h(0). □ 

How do mutually unbiased bases fit into this scheme? The answer is 
that they almost perfectly do, except for the fact that the underlying 
group is not Abelian. Indeed, let G = Udxd the group of d x d unitary 
matrices, and let H C Udxd denote the set of complex Hadamard ma- 
trices (rescaled by the factor 1/yd)) in Udxd- Let the 'forbidden' set 
A be the complement of H. Of course, the group operations + and — 
in the Delsarte scheme are now replaced by matrix multiplication and 
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inverse. Also, the role of zero element is taken by the identity matrix. 
Then, the maximal number of mutually unbiased bases in C d is exactly 
the maximal cardinality of a set {Uq, Ui, . . .U m } C G such that all 
'differences' U*Uk (0 < j,k < m) lie in the prescribed subset A c U {/}. 

Unfortunately, we do not know how to generalize Delsarte's method 
to the case of non-commutative groups, in particular to G — U dX d- Nev- 
ertheless, we can still use Delsarte's scheme if we rephrase the problem 
appropriately, as follows. 

Assume that a family Hi, . . . H m of m mutually unbiased complex 
Hadamard matrices exists. Then all entries of all matrices are of mod- 
ulus 1, and the columns (and thus the rows) within each matrix are 
complex orthogonal, and we have the unbiasedness condition: for any 
two columns u, v coming from different matrices we have | (u, v) | = \fd. 
(Recall that we have re- normalized the matrices by a factor of y/d.) 

After multiplying rows and columns by appropriate scalars if nec- 
essary, we can assume that all coordinates of the first row and col- 
umn of Hi are l's, and all coordinates of the first row of all other 
matrices are l's (i.e. we assume that all appearing columns have 
first coordinate 1, and the first column in Hi consists of l's. This 
is standard and trivial normalization.) All the other coordinates in 
the matrices are complex numbers of modulus 1, i.e. they are of the 
form e 2wip with p e [-1/2,1/2). Therefore, we can associate to each 
column vector (1, e 27rip1 ,. . . , e 2 ™^- 1 ) the vector (0, pi,..., p d -i) € T d , 
the real d- dimensional torus, T d = [—1/2, l/2) d . Also, note that the 
first coordinate always automatically becomes 0, because each column 
starts with coordinate 1. Therefore we make the more useful asso- 
ciation that a column c = (l,e 2mpi , . . . ,e 2mpd - 1 ) is represented by 
u = (pi, . . . Pd-i) € T^ -1 , the d — 1-dimensional torus. There are 
altogether md column vectors in the Hadamard matrices Hi, . . . H m , 
and we will denote the associated vectors in T^ 1 by Ui, . . . u md (we 
will see that in this approach it is not really relevant to indicate which 
vector comes from which basis. But let us agree for convenience that 
ui = (0, . . . , 0), corresponding to the first column of Hi.) 

Two columns Ci = (1, e 2nipl , . . . , e 2 ™^- 1 ) and c 2 = (1, e 2nipl , e 2 ^- 1 ) 
are orthogonal if and only if 1 + Yl'jZi e 2m ^~ Pj+,J '^ = 0, and they are 

unbiased if and only if [1 + e 27r ^ _ft+M ^| = Vd. Therefore it is 

natural to introduce the following definitions. 

Definition 2.1. Let ORT d denote the set of vectors (ai, . . . a d -i) € 
T d_1 ; in the d — 1-dimensional torus, such that 1 + Yl'jZi e 2m(X] = 0. 
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Also, let UB d denote the set of vectors (a±, . . . a d -i) G T d l , such 
that |1 + Xlj=i e 27rlQ!i | = \fd. Let us also define the 'forbidden' set 
A d = (ORT d UUB d y. 

We conclude that the vectors u 1; . . . u md satisfy that the difference of 
any two of them (the difference being taken mod 1 in each coordinate, 
i.e. we take the difference in the group T d_1 ) lies in A^U{0}. Therefore, 
we have arrived exactly to the scheme of Lemma 12. f I 

As a preliminary remark we note that the dual group of G = T d_1 is 
G = Z d_1 . And the action of a character 7 G Z d_1 on a point x G T d_1 
is given as 7(2;) = e 27r ^ 7 '^. In particular, 7 = is the trivial character 
(constant 1). The Fourier transform of a function / : G — > C is a 
function / : G — > C given as f(j) = f G f(x)j(x)dx. 

Let us see whether we can find a good 'witness' function in this 
situation. At first sight things do not look promising because we have 
no understanding of the geometry of the sets ORT d and U B d inside the 
torus TP* -1 . However, it turns out such geometric understanding is not 
necessarily required and we easily prove the following generalization of 
Theorem 11.11 

Theorem 2.2. Let A be an orthonormal basis in C d , and let B = 
{ci, . . . c r } consist of unit vectors which are all unbiased to A. Assume 
that for all 1 < j 7^ k < r the vectors Cj and c& are either orthogonal 
or unbiased to each other, i.e. either (cj, c&) = or |(c 3 -,Cfc)| = l/y/d. 
Then r < d 2 . 

Proof. As we saw in the discussion above, the vectors u 1; . . . u r G T^ 1 
(associated to yfdci, ■ ■ ■ \fdc r ) satisfy Uj — u k G A c d U {0} for all 1 < 
j, k < r. Therefore Lemma 12.11 can be applied. 

Define the 'witness' function h : T^ 1 — y R as follows: 

d-l 

-£< 

It is trivial to check that h satisfies all requirements. Indeed, h is an 
even function which vanishes on ORT d U UB d . The Fourier coefficients 
of h are simply the coefficients of the exponential terms after expand- 
ing the brackets, and these are clearly nonnegative. Also h(0) = 1 
because h(0) is the integral of h, which is just the constant term. Also, 
h(0, . . . 0) = d 2 , so that we conclude from Lemma [2~T1 that \B\ < d 2 . □ 



h(x 1 , . . . Xd-x) 



(d-l)d 



\i + j2 e2mX] \ 2 - d 

3=1 
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Remark 2.3. As shown by Theorem 11.21 the result of Theorem 12.21 is 
sharp if d is a prime-power. However, if d is not a prime-power, then 
it could be possible to find a better witness function than above. The 
function h above uses simply the definition of the sets ORT d and UB d . 
In principle, it could be possible to find some structural properties of 
these sets in dimension 6 (or any other composite dimension), in order 
to construct a better witness function and get a sharper bound on r. 
Any upper bound r < d 2 would mean that a complete set of MUBs 
does not exist in dimension d. We have not been able to make such 
improvement for any d so far. □ 

Another observation is that if r = d 2 in Theorem 12.21 then both 
estimates (j2J), (jSJ) must hold with equality. On the one hand, it is trivial 
that (j3J) automatically becomes an equality for the h above (because h 
is zero on ORT d and UB d ). On the other hand, inequality (J2J becomes 
an equality if only if \B{pf)\ 2 h{pf) = for all 7^0. These are non- 
trivial conditions and we obtain the following corollary, which is a 
generalization of Theorem 8 in [3] . 

Corollary 2.4. Let A be an orthonormal basis in C d , and let B = 
{ci, . . . Cfp} consist of unit vectors which are all unbiased to A. Assume 
that for all 1 < j ^ k < d 2 the vectors Cj and are either orthogonal 
or unbiased to each other. Write B as a d x d 2 matrix, the columns of 
which are the vectors Cj, j = 1, . . . d 2 . Let ri, . . . r^ denote the rows of 
the matrix B, and let r^ = r^/r^ denote the coordinate-wise quotient 
of the rows. Then the vectors r^ (1 < j > j£ k < d) are orthogonal to 
each other in <C d2 , and they are all orthogonal to the vector (1, 1, ... 1) G 

Proof. This is a direct consequence of the proof of Lemma [2TTT Indeed, 
for (T5]) to be an equality -6(7) must be zero whenever h(j) 7^ and 
7 7^ 0. Looking at the definition of h, ^(7) 7^ happens exactly when 
7(^1, • • • Xd-\) = e 2m ^ X3 ~ Xk)+ ( Xq ~ Xs ' > \ for any quadruple < j, k, q, s < 
d — 1, where we use the convention that xq = 0. Using the definition 
of B(j) we obtain that -6(7) = means exactly that r^y* and r s / q are 
orthogonal to each other. 

□ 

Remark 2.5. There are altogether d{d— 1) vectors of the form r^, and 
with the addition of (1, 1, ... 1) we get a system of d(d— 1)+1 orthogonal 
vectors in C d2 . It is not at all obvious whether in each dimension d 2 
there exists a set of vectors R = {r±, . . . r^}, with all coordinates having 
absolute value 1, such that they satisfy these orthogonality constraints. 
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Let us remark that if modulo d 2 a <i-element Sidon-set exists (i.e. a set 
where the non-zero differences are all distinct), then the corresponding 
rows of the Fourier matrix F d 2 form an appropriate system R. For 
example, for d = 6 the set {0, 1, 3, 8, 23, 27} is a Sidon-set modulo 36, 
so the rows R — {f , fi, f3, fs, f23, ^27} of the Fourier matrix F 36 satisfy all 
orthogonality constraints (but R is not a concatenation of orthonormal 
bases, so it does not form a complete set of MUBs). One could try to 
find a composite d such that a <i-element Sidon set modulo d 2 does not 
exist, and then show that an appropriate set of vectors R also cannot 
exist. This would prove that a complete set of MUBs does not exist in 
dimension d. □ 

3. Linear duality and pseudo-MUBs 

We can view the set B in Lemma [2. II as a 0-1- valued function on G. 
Also, observe that B does not directly enter the proof, but instead the 
function |B(7)| 2 = B — -8(7) is essential. For any y G T rf_1 let f(y) 
denote the number of ways of writing y as a difference of two elements 
of B. Then, for any 7 G G = 7L d ~ x we have | (t) 1 2 — J2j k 7(fy — ^fc) = 
Yl y f{y)e 2wi ^' y ' v '- Therefore, / has the following essential properties. 

• the finite exponential sum ^ y f{y)e 2m ^ l,y ' is nonnegative for all 
7 G Z d -\ and the exponents y G A c U {0} = ORT d U UB d U {0}. We 
can view it as the Fourier transform of the function / : T d_1 — > M + . 

• the coefficients f(y) are nonnegative integers. 

• the sum of the coefficients f(y) = \B\ 2 - 

• /(0) = \B\. 

Given any such function / we can repeat the proof of Lemma 12.11 
with the 'witness' function h defined in (j3J), and conclude that 

' /(o) - ft(0) 

This motivates the following definition: 

Definition 3.1. We will call a function f : T^ 1 — > R + a complete- 
pseudo-MUB-system in dimension d (or pseudo-MUB-d in short) if it 
satisfies the following conditions: 

• f is nonnegative, the support of f is finite and is contained in 
ORT d U UB d U {0} 

• the finite exponential sum f(j) = J2 y f(y)e 2m( -' y ' y ' ) is nonnegative 
for all 7 G Z^ 1 . 
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• the sum of the coefficients Ylyfiv) = ^ ■ 
•f(0) = d 2 . 

Notice that / is not required to be integer valued. As discussed 
above, a complete set of d + 1 MUBs always gives rise to a complete- 
pseudo-MUB-system. Indeed, let B C T d_1 denote the d 2 columns 
of the corresponding d mutually unbiased Hadamards, and let f(y) = 
(B—B)(y), meaning the number of ways y can be written as a difference 
bj — bk- Then / is a complete-pseudo-MUB-system. The converse is 
not necessarily true: a complete-pseudo-MUB-system does not directly 
imply the existence of a complete set of MUBs. 

In any dimension d, if we find a pseudo-MUB-d function / then it 
could serve as a 'dual- witness' testifying that our function h in equation 
(jlj) is best possible, and it would mean that the Delsarte method alone 
cannot prove the non-existence of d + 1 MUBs in dimension d. We 
emphasize that it would not mean that a complete system of d + 1 
MUBs exists. It would only mean that a complete-pseudo-MUB-system 
exists. 

We remark that there is a linear duality here. Either a pseudo-MUB- 
d exists, or a better witness function h : T d ~ 1 — > R exists, proving that 
r < d 2 in Theorem 12.21 There is no third option! In the latter case 
we could conclude that no complete system of d + 1 MUBs exists in 
dimension d. In the former case we would have an interesting pseudo- 
MUB-<i in our hand, which could possibly lead later to the discovery 
of a proper complete set of MUBs. 

Let us examine the situation in dimension d = 6. 

One natural idea is to fix some m, and look for a pseudo-MUB-6 
function such that its support contains vectors only whose coordinates 
are mth roots of unity. The reason is that all known complete sets 
of MUBs consist of such vectors. It is also convenient because such 
vectors belonging to ORT 6 and UB Q can easily be listed by a computer 
code. Furthermore, the restriction /(7) > needs to be checked only 
as 7 ranges over the cube [0,m — l] 5 , due to periodicity. Finally, we 
fix /(0) = 1 (which is a somewhat more convenient normalization than 
/(0) = d 2 in the definition), and maximize M = Ylyfiv) by linear 
programming (a pseudo-MUB-6 would have the value 36 here). We 
have tried this and the results are the following: 

• m = 12, M = 17.5 

• m = 8, M = 21.6 
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• m= 16, M = 21.6 

Larger values of m are out of our computational power. As one 
can see, these results are inconclusive. We could not find a pseudo- 
MUB-6 but we could not find a better 'witness' function h(x) either in 
dimension 6. By linear duality one of them must exit, and it would be 
interesting to see which one. 
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